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ABSTRACT 
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We analyze the evolution of scalar and gauge fields during first order phase 
transitions and show how the Kibble mechanism-'^) for the formation of topologi- 
cal defects emerges from the underlying dynamics, paying particular attention to 
problems posed by gauge invariance when a local symmetry is spontaneously bro- 
ken. We discuss also the apphcation of the mechanism to semilocal defects and 
electroweak strings. 
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1. Introduction 



Topological defects arising in theories with a spontaneously broken global 
or local symmetry play an important role in various branches of physics from 
cosmology^'^) to condensed matter physics'^ 

In cosmology, topological defects are expected to form during phase transitions 
in the early Universe. It is important to know the number density of defects at 
the time of formation, as well as other statistics such as correlation functions. 
Depending on the initial number density, many physics models admitting local 
monopoles or domain walls are ruled out®''. For defect-induced baryogenesis^'' it 
is important to know the initial separation of defects, and for structure formation 
scenarios using topological defects - cosmic strings^^^ , global monopoles^^^ or global 
textures-*^^) - it is also useful to understand the initial distribution of defects. 

In 1976, Kibble^) suggested a simple mechanism for the formation of defects in 
cosmological phase transitions, where the order parameter is a scalar quantum field 
(p (which may be elementary or composite). It has two key ingredients: randomness 
of the phases of the scalar field on length scales larger than some correlation 
length ^, and the geodesic rule for interpolating the values of on curves connecting 
points in different correlation volumes. The first assumption states that takes 
on random values in its vacuum manifold Ai at points separated by a distance 
greater than ^, the second one states that along a curve in space connecting two 
such points, the field (f) traces out a geodesic path on the vacuum manifold. 

For global defects, the Kibble mechanism has been verified in liquid crystal 
experiments^^^ and in numerical simulations in field theory^^^. For local defects, 
convincing experiments are lacking (see, however, Refs. 15 for some numerical 
results for local vortex formation) . 

Recently, Rudaz and Srivastava^^-* have argued that in local field theories, the 
geodesic rule is not justified and that the rate of defect formation might be much 
smaller than what would be obtained by naively applying the Kibble mechanism. 
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The core of the objection is that the Kibble mechanism is formulated in a gauge- 
dependent way, and that a geodesic curve in the vacuum manifold of (j) does not 
necessarily minimize the energy density and hence should not play a distinguished 
role. These objections highlight the need to reconsider the Kibble mechanism for 
defect formation in theories with local symmetry breaking. 

In this paper we reconsider the Kibble mechanism for global and local defect 
formation at first order phase transitions. We analyze the equations of motion for 
scalar and gauge fields, and demonstrate the validity of the geodesic rule for 0. To 
be specific, we consider vortex (in 2 + 1 dimensions) or cosmic string (in 3 + 1 di- 
mensions) formation in a model in which a [/(I) global or local symmetry is broken. 
However, our methods also work for more complicated models with 7ri(A1) ^ 1, for 
other types of topological defects, and for semilocal^^^ and electroweak^^) strings. 
We do not here consider second order transitions because the analysis is qualita- 
tively different: the semiclassical methods we use ignore thermal fluctuations. 

In Section 2, we discuss theories with a first order phase transition which 
procedes via bubble nucleation. We discuss how far the analysis can be applied 
to semilocal defects and electroweak strings in Section 3, and Section 4 contains a 
brief summary of results. Throughout the paper, units in which ks — Ti — c — 1 
are employed. Greek indices run from to 3 and our space-time signature is 

(+,-,-,-)■ 

2. First order phase transitions 

To be specific, we shall consider an Abehan Higgs model with a complex scalar 
field and a [/(I) gauge connection A^. The Lagrangian of the system is 

C = {D^(t^)^D^(t^ - V(0) - i F^^Fi'- (2.1) 

where 

D^^df,- ieA^ (2.2) 
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is the covariant derivative (e is the gauge couphng constant), 



F^^ = d^Ay - d^A^ (2.3) 

is the field strength, and V^(0) is a symmetry breaking potential, whose zero tem- 
perature form is A(|0p — rf')'^ /2. We assume that T^(0) depends only on |0|. In 
this case, the Lagrangian (2.1) has a local C/(l) symmetry. If |0| = 77(T) 7^ is the 
absolute minimum of y(0), then this symmetry is broken at low temperatures. We 
shall assume that = is the minimum of the high temperature effective potential 
Vt(0). At a critical temperature Tg, = ceases to be the absolute minimum of 
Vt{^). 

Considering the Abelian Higgs model in the context of an expanding Universe, 
we conclude that as the temperature drops below T^i the system undergoes a 
symmetry breaking phase transition. The order of the transition depends on the 
functional form of ((/?). It is generally thought that the transition is first order 
for (3 = A/e2 ^ 0, and may be second order for (5 ^ 00}^^ Iti the former case, the 
one-loop temperature dependent effective potential has the form^^^ 

Vm = \x{Tm^ - rff + -/T''\^\^ - ^e'm\ (2.4) 

In a first order phase transition, = remains a metastable fixed point below 
T — Tc- The transition to a state with |0| = 77 procedes by bubble nucleation.^-'^^ 
There is a finite probability per unit volume per unit time dP/dVdt that a bubble 
with = r] will nucleate in a surrounding sea of "false vacuum" = 0. This 

probability is given by an expression of the form^-*^'^^) 

dP 

Aexp-SE[<j>,A] (2.5) 



dVdt 



where and A are extrema of the Euclidean action Se- This field configuration 
represents a tunnelling process. At high temperatures (high compared with the 
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scalar and vector masses nis and my) the action is effectively three dimensional, 
and the tunnelling solution is spherically symmetric about a point. ^^-^ If the energy 
difference between the minima at (j) — and \4>\ = r]{T) is small compared with 
the height of the barrier between them, then the solution is well approximated 
by a thin walled spherical bubble of true vacuum inside the false one.^^^ The 
width of the bubble wall is approximately m~^, and the radius of the bubble is 
of order Vb/mgAT/, where Vj, is the barrier height and AT/ is the difference is 
free energy between the two phases. The bubble will then expand with a speed 
which depends on the interaction between the wall and the rest of the hot matter 
in the universe. We shall assume that the couplings are small enough that we 
may take the limiting velocity to be 1. The expansion of the bubble is fueled 
by the conversion of potential energy density V{Q) to wall kinetic and gradient 
energy. The phase transition is completed when neighboring bubbles collide and 
the fraction of space with |0| = 77 approaches unity. The correlation length ^ for 
this transition is the mean separation of bubbles. Implicit in our analysis are some 
assumptions about time scales: we are taking the expansion rate to be much longer 
than the nucleation rate per Hubble volume, so that we are justified in taking the 
spacetime to be Minkowsi. The expansion of the universe then serves to reduce 
the temperture as a known function of time. 

According to the Kibble mechanism, there is a fixed probability p of the order 1 
that a defect will form in any correlation volume The exact value of p depends 
on the type of defect, i.e., on the topology of the vacuum manifold M. (see e.g., 
Ref. 23 for a recent calculation of p for various models). 

Let us illustrate the Kibble mechanism for our toy model, the Abehan Higgs 
model. The phase a of is assumed to take random values in different bubbles. 
After two bubbles meet, then if we follow a line in space connecting the centers 
of the two bubbles, a is assumed to interpolate more or less smoothly between its 
values in the two bubbles. The second part of Kibble's argument states that a{x) 
will follow a geodesic in M. in order to minimize the potential energy. 
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A vortex can form when three bubbles colhde - as ilhistrated in Fig. 1. If 
ai, i = 1,2,3 are the phases of in the three bubbles, then a vortex will form if 
the sum of the phase differences 0:2— and 0:3— 02 exceeds vr. In this case, a{x) will 
run from to 27r as we go along the circle 7, i.e., the field configuration has winding 
number 1. There is then a topological obstruction to the scalar field reaching the 
vacuum manifold everywhere in the region bounded by the lines connecting the 
centres of the bubbles. The field must vanish somewhere, and this is where the 
vortex forms. 

As discussed in Ref. 24, a bubble collision is a quite violent event. In order 
to justify the Kibble mechanism (both for global and local symmetry breaking), 
we must follow the evolution of amplitude and phase of (f) and demonstrate that 
the geodesic rule is valid. In order to do this, we will write down the dynamical 
equations which follow from (2.1). We first consider a global theory and discuss 
the evolution of amplitude and phase of 0. Then, we make the transition to a local 
theory and show that the gauge fields do not strongly perturb the evolution of 0. 
This analysis is done in a particular gauge. However, the final result - the winding 
number - is gauge invariant. 

In Lorentz gauge, i.e., with the choice 

d^Ai" = (2.6) 
the variational equations which follow from (2.1) are 

(9^9^ - 2ieA^d'' - e^A^A^) <j> + 2 ^<l> = (2.7) 

and 

d/^d^'Ay - 2e^Ay\(t)\^ = -ei(l)*dy(f) . (2.8) 
It is convenient to separate Eq. (2.7) into equations for the amplitude p and phase 
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a oi (J). Inserting 



(2.9) 



into (5), we obtain 

dV 

d'^p-{da-eA)^p-e^A^p + 2-^p = (2.10) 

and 

d^a + 2{d''a-eA'')d^p^^Q. (2.11) 

The collision of two bubbles in the Abehan Higgs model was studied numeri- 
cally in Ref. 24. We now demonstrate that we can reproduce the essential features 
of the collision process using the above equations. 

We first consider a theory with global U{1) symmetry and choose axes such that 
the centers of the bubbles lie on the x axis. We can additionally boost in the (y, z) 
plane to a frame in which the bubbles nucleate simultaneously^^^ , and translate in 
the X direction so that they collide at x = 0. Provided the bubbles nucleate far 
from each other (far meaning much greater than the wall width so that the field 
is exponentially close to zero between the bubbles) the solution representing two 
bubbles nucleating at xi and X2 can be approximated by a sum ansatz 

0(f , 0) = e'^'fix - f i) + e^°V(^ - X2) (2.12) 

where / is the modulus of the single bubble field. The action for this configuration 
is (exponentially) independent of ai and a2'- hence the phases of the two fields are 
well approximated by independent random variables. Without loss of generality 
we take the phase in one bubble to be a = and in the other bubble a — ac- Note 
that the phase a is not defined for points outside the bubbles. 
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When the bubbles meet (we take this to occur at time t = 0), the phase a{x) 
is approximately a step function 



a{x,t = 0) = ac9{x) . (2.13) 

In fact, the step will be smoothed on a scale by the finite thickness of the 
bubble walls. We shall work in the planar approximation dy(j) = dz4> = 0- This is 
reasonable if the bubble radius is much greater than the wall thickness. In fact, 
the symmetry of the two bubble collision reduces the problem to a two dimensional 
one anyway^^) , but for simplicity we choose not to exploit the coordinates in which 
this is manifest. Inserting (2.13) as initial condition into the phase equation (Eq. 
(2.11)), we see that phase wave fronts emerge which travel in ±x direction with 
the speed of light (see Fig. 2): 

a{x, t>0)^ \ac [9{x + i) + Q{x - t)\ . (2.14) 

Eq. (2.14) gives the solution of (2.11) because the phase waves are propagating 
inside the bubbles where ^^J,p — 0. 

The phase waves which arise for etc 7^ carry away some of the kinetic energy of 
the walls, but not all. The rest of the energy goes into bubble walls. The region of 
false vacuum = does not disappear at t = 0. Rather, the walls which separate 
the region with \<^\ = ri from the false vacuum scatter and start to reexpand^*^) . 
What is happening is that the modulus of the field has overshot p = y], and rolled 
back to p = 0. We denote the wall positions by ±X{t). However, as X{t) increases, 
the potential energy of the field configuration increases, thus creating a force which 
tends to restore X to 0. An approximate equation for X{t) can be obtained from 
elementary physics considerations: 

dW 

.X = -— , (2.15) 

where a is the mass per unit area of the wall, and W{X) is the potential energy 
(per unit area) of the false vacuum region x e [—X, X\ . Obviously, W — XV {Q) 



9 



and hence (2.13) becomes 

X = -T/(o)/(7 (2.16) 

Thus, there is a constant restoring force which causes the wall to recollapse on 
a time scale r = 2X{0)a/V{0) . There is a series of bounces, each losing energy 
through propagating oscillations in p, and thus leading to reduced ^(0) and period 
T (see Fig. 2). 

An alternate way to derive Eq. (2.15) is to insert the ansatz 

p(x,t)^r]e{x-X(t)) (2.17) 

into the equation for the modulus p (Eq. (2.10)), and integrate the resulting 
distributional equation over x. 

To summarize, we have verified that after a bubble collision, the phase a{x) 
along a path 7 connecting the two bubble centers interpolates between its original 
values in the two bubbles. The interpolation happens in two jumps of width 
associated with excitations of the Nambu-Goldstone mode, spreading from the 
collision site at the speed of hght. Thus, the second key ingredient of the Kibble 
mechanism, the "geodesic rule," has been estabhshed for global defects forming in 
a first order phase transition. We stress that the geodesic rule follows from the 
equations of motion, not from minimizing the energy, as assumed in Ref. 16. 

We now extend the analysis to theories with a local U{1) symmetry. We again 
study the collision of two bubbles and establish the applicability of the geodesic 
rule. In the bubble nucleation instanton, the magnetic field associated with 
vanishes. Thus at the time of nucleation in we have, even in the Lorentz gauge, 
some freedom left in specifying A^. In general, 

A{x, Q = VA(f , in), = A(f , in), (2.18) 

with A{x, in) an arbitrary function of position. However, it is clearly simplest 
to choose A(x, in) = = A(x, in), which amounts to a complete specification of 
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the gauge. During the bubble wall collision, a nonvanishing gauge connection is 
generated through the coupling to the phase difference in the scalar field across the 
bubble wall, which in turn feeds back into the evolution of the scalar field modulus 
p via (2.10). We will analyze firstly assuming A/e^ S> 1. We assume that we may 
neglect eA in comparison to da ~ ^^\r], so that Eq. (2.8) becomes 

(a| - dl)Ay = 2e|0|2 {dya - eAy) ~ 2erj^dya . (2.19) 

If we take for a Eq. (2.14) then we can easily solve for the gauge field. This is 
most easily expressed in light cone coordinates = t ± x. After one integration 
we find {t > 0) 

d-A+ = leri'^ac9{x+), 

(2.20) 

d+A- = -ler)'^ac9{x ), 
which corresponds to an electric field in the x direction, of 

Foi = 2(d-A+ + d+A_) = ^'n'^adOix + t)- e{x - t)]. (2.21) 

This field exists only between the phase waves, which are carrying off equal and 
opposite charges (per unit area) ad 2. Behind the phase wave the Higgs vacuum 
screens the charges, and the gauge field is exponentially damped beyond lightcone 
distances Ax^ ~ [er])-^. Thus we conclude that the gauge field is bounded in 
modulus by ~ i]. Our approximation is self-consistent, for er] <^ v^rj. At the 
phase wave itself the gauge field vanishes, so it should be a good approximation to 
ignore its effect on the propagation of the wave (this feature was originally found 
by Hawking et al in their numerical simulations^^-*). 

This is perhaps not the correct limit to use if we are assuming a first order 
phase transition, where A/e^ is supposed to be small. In this limit we cannot 
approximate the phase wave by a step function. Instead we take {t > 0) 

«W = iMl + iH/Kx+) + ^M/K.-)l (2.22) 

where W is the phase wave profile, interpolating between -1 and 1 as its argument 
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changes from —oo to +00. Then we find that Eq. (2.8) becomes 



(2.23) 



In a gradient expansion in powers of mg/my, the first term is A± ~ d±a/e — 
:^acmsW' {rrisX^) / A. Thus the current vanishes to this order, and the equations of 
motion for a and p are affected only at higher order in mg/my. 

Physically, what is happening is that in the limit mg ^ my, the wall collision 
contains enough high frequency modes to create longitudinal gauge bosons. In the 
opposite limit, the gauge bosons are too massive to appear from scalar bosons of 
freqency ~ mg: instead, the gauge field tracks the phase to ensure that the current 
vanishes. 

In either case, the key point is that the interpolation of the phase between the 
bubbles is not affected by the presence of Ay. Hence, by chosing a particular gauge, 
we establish the geodesic rule for local theories. In a three bubble collision, the 
geodesic rule can be applied between each pair of bubbles provided the collision 
happens before the advancing wall of the other bubble reaches the collision point. 
We can in fact boost along the normal to the plane containing the centres of the 
three bubbles so that they nucleate simultaneously, so this is a constraint on the 
positions of the centres: they must form an acute triangle. The winding number 
around a closed curve connecting these points, along which the scalar field vanishes 
nowhere, is a gauge invariant quantity. This is given by 



We can compute it in our chosen gauge and be confident that the result is gauge 
invariant. 

These results imply that in first order phase transitions, there is a finite proba- 
bility that a field configuration of nontrivial winding (and hence a vortex) emerges 




(2.24) 
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during the collision of three bubbles. This is true for both local and global the- 
ories. If ai, ai and ^3 are the phases of in the three bubbles (without loss of 
generality a\ < a2 < as) and if ai + tt < as < 0:2 + tt, then by the geodesic 
rule, after bubble collision, a{x) will smoothly go from to 27r along 7, yielding 
a configuration with winding number 1. As we have seen, the phases in the three 
bubbles are random, since the nucleation probability is independent of their values, 
and so we get winding number 1 (with probability 1/4.^^)) 

3. Extensions 

Our analysis has been based on dynamical rather than on topological con- 
siderations. Hence, the arguments may carry over to the case of defects such as 
semilocal strings^ ''^^ or electroweak strings^^^, whose stability or otherwise is dy- 
namical in origin. 

Key to our analysis in previous chapters was to establish the validity of the 
geodesic rule, i.e., of the statement that after completion of the phase transition, 
the phases of the Higgs fields along a line in space connecting the centers of two 
initial bubbles will interpolate smoothly between their initial values in the bub- 
bles, thus forming section of geodesies on the vacuum manifold M.. The question 
of determining the probability of defect formation then reduces to the statistics 
problem of finding the probability that a closed geodesic will have nonvanishing 
winding number. 

We now generalize this approach to semilocal and electroweak strings. We first 
verify the geodesic rule, thus reducing the problem to a probability problem. How- 
ever, the probability calculation will in general be much harder than for topological 
defects. 

Semilocal strings^^-* arise in models with a large global symmetry group G of 
which only subgroup G\ is gauged. If this gauge group and its unbroken subgroup 
B.\ obey the usual topological condition t:i{G\/H\) ^ 0, then stable string solutions 
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exist only if the scalar mass is small enough relative to the vector mass. 

26,27) Tj^g 

simplest example is the model of Ref. 17: a U(2) symmetry acts on a complex 
scalar doublet but only the U(l) generated by the identity matrix is gauged. 
When $ gains an expectation value the remaining symmetry is a global U(l). 

The vacuum manifold of the theory is S^, defined by = 77^, and it is fibred 
by the action of the local U(l) into a bundle which is locally S"^ x - the Hopf 
bundle. In the string sohition, the scalar field wraps around one of these fibres 
outside the string, and vanishes at the origin. In fact, it is simply an embedding of 
the Nielsen-Olesen^^) vortex in the full theory. It is only stable if the scalar mass 
is less than the vector mass. 

The equations of motion of this theory are essentially the same as those of 
Section 2. Both components 0i and 02 of $ satisfy Eq. (2.5): 

dV 

d^(t>i - 2ieA^9'^0i - e'^A^(t>i + 2 ^ 0^ = i = 1, 2 , (3.1) 

and the gauge field satisfies the analog of Eq. (2.8): 

d^Ay - 2e^Ay{(t)l + 0i) = -ie{(t)ldy 0i + 0^9^ 02) . (3.2) 

In particular, the gauge connection does not mix 0i and 02- The same argu- 
ments as in Section 2 imply that the geodesic rule applies, in our Lorentz gauge 
with A^(x, t-a) — 0. 

However, this is a geodesic rule in the full vacuum manifold 5"^, in which there 
is no topological obstruction to the scalar field reaching |$| = everywhere in the 
triangular region formed by a three-bubble collision. This is related to the fact 
that the analogue of (2.24), or 

is not necessarily an integer. The probability of forming a semilocal string is in 
fact a complicated dynamical question. It seems likely however, that the 'closer' 
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the field in the three bubbles lies to the same U(l) orbit, the more likely is the 
formation of a string. We can in fact give this closeness a precise geometrical 
meaning. 

Let us choose coordinates (cr, i>,x) on such that in its vacuum manifold 



$ = ?7 



cosx/2 e' 
sinx/2 e«(^+V')/2y 



(3.4) 



Then a parameterizes the fibres, and (x, '4') are polar coordinates on the base 
space S'^. This can be seen by projecting onto a unit 3- vector 



(p"' — $t(ja$/|$|2 _ (sin X COS ■0, sin X sin ^0, COS x)- 



(3.5) 



The semilocal string has (x, ip) constant around it, with a changing by An. Thus 
the measure of closeness of the field values in the three bubbles (crj, Xi, V'i) to the 
string configuration is the area of the spherical triangle defined by {xiii^i)- The 
smaller this area, the closer the field is to a pure phase change around the curve 
joining the centres of the bubbles (that is, the closer is to 1), and the more 
likely a string is to form. Unfortunately, without dynamical simulations, we cannot 
say how the string formation probability depends on the area A. The only piece 
of information we can extract geometrically is the probability P{A) that a random 
spherical triangle has area less than or equal to A, which is^^^ 



P{A)- 
For small A, 



A 
2" 



A A 



A 



— sm — cos TT sm — 

2 2 2 2 



^ ' 47r 



i(7r-^)(37r-A)tan| 



/ 2^ 
/TTCOS — . 

' 2 

(3.6) 



(3.7) 



For fields $ with d components we would expect the probability to go as V , where 
y is a 2d—l dimensional volume, although the coefficient is a more difficult exercise 
in geometric probability. 
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For electroweak strings the geometrical considerations are identical: our two 
component semilocal model is just the bosonic sector of the electroweeak theory in 
the limit that the weak mixing angle approaches 7r/2. However, we can expect the 
fields to evolve towards a string solution only if there is a locally stable solution 
towards which to evolve. In the Standard Model the string is unstable^^^ , so we do 
not expect to see them formed. However, metastable strings exist in theories with 
more complex Higgs sector, and in those we can expect some string formation 
in first order phase transitions, although we are unable to estimate the probability. 

4. Conclusions 

We have studied the formation of topological defects in first order phase tran- 
sitions. We showed that both for local and for global defects, the assumptions on 
which the Kibble mechanism^) is based can be established by using the equations 
of motion. 

The results of the analysis is that there is a probability p that a defect will form 
per correlation volume of the field. For topological defects, p is of the order 1^^), for 
nontopological defects such as semilocal strings and electroweak strings, p depends 
on as yet ill-understood dynamics, but we guess it to be typically much smaller than 
1. However, in no circumstance is the formation probability Boltzmann suppressed. 

Our technique is based on first studying the dynamics of the scalar fields in 
the absence of gauge fields (a gauge dependent analysis), establishing the existence 
of winding number in the final state (a gauge independent conclusion), and then 
bounding the effects of the gauge fields, for convenience working in Lorentz gauge 
with the additional choice ^4^ = at the time of bubble nucleation. 

We have not included thermal fluctuations into our consideration. We are 
therefore implicitly assuming that the root mean square amplitude of the thermal 
fluctuations on the length scales we discuss is much less than the magnitude of the 
scalar held r]{T), which is in fact not unreasonable.^^^ 
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5. Figure captions 



Figure 1 

Three bubbles of the broken symmetry phase (p = rj) colhding. If the phase 
change of the scalar field around the loop 7 is ±2%, a string (or antistring) is 
formed. If the phases ccj are ordered, then then the requirement for a string is 

CKl + TT < 0:3 < Q;2 + TT. 

Figure 2 

Space-time diagram of two bubbles with different phases colliding. After nu- 
cleation at time tn the bubble walls collide at t = 0, when they are travelling at 
approximately c. Phase waves continue out from the collision site, spreading a 
region whose phase is halfway between the phases of the bubbles. The walls pass 
through each other, but eventually turn round and recoUide. This may happen 
several times. 
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